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The recent experimental conductance measurements taken on magnetic impurities on metallic 
surfaces, using scanning tunneling microscopy technique and suggesting occurrence of inelastic scat- 
tering processes, are theoretically addressed. We argue that the observed conductance signatures 
are caused by transitions between the spin states which have opened due to e.g. exchange coupling 
between the local spins and the tunneling electrons, and are directly interpretable in terms of inelas- 
tic transitions energies. Feasible measurements using spin-polarized scanning tunneling microscopy 
that would enable new information about the excitation spectrum of the local spins are discussed. 



Inelastic scattering processes play a crucial role in 
studies of excitation spectrums of e.g. single atoms and 
clusters of atoms, molecules, quantum dots, carbon nan- 
otubes, and graphene. Signatures of inelastic scattering 
are often displayed in conductance measurements as steps 
that cannot be ascribed to the resonant elastic tunneling, 
something which led to the discovery of molecular vibra- 
tions in metal-insulator-metal junctions exposed to pro- 
pionic acid [1] . Inelastic features arising from vibrations 
have since been studied in e.g. molecular structures [2J, 
Josephson junctions , and surfaces |3] . 

Recently, inelastic scattering have been observed 
in many studies of magnetic structures. Magnetic 
anisotropy of magnetic atoms and clusters have, for in- 
stance, been studied through spin inelastic scattering 
[3J m IZl IHl in]- Signatures of spin inelastic scattering 
processes are, furthermore, frequently appearing in ex- 
citation spectrums taken on bulk materials, e.g. Fe, and 
thin films, e.g. Co [10], and individual magnetic adatoms, 
e.g. Fe and Cr [11], Co [El IE], and molecular magnets 
|14| . In contrast, there are only a few theoretical studies 
in which spin inelastic scattering processes are included 
in a non-equilibrium description |15l 1161 117j . 

Here, we address the theoretical description of inelas- 
tic signatures in the (differential) conductance through 
a tunneling barrier in which a local spin is embedded, 
and we include the possibility of magnetic source and 
drain. The set-up is, thus, pertinent for scanning tun- 
neling microscope with or without spin-polarized tip. As 
we discuss below, the conductance can be naturally di- 
vided into three terms of which the last is sensitive to 
the spin fluctuations in the local spin structure. Those 
fluctuations include transitions between the ground state 
and excited states which are forbidden in the absence of 
interactions between the local spins and the surrounding 
electrons in e.g. tunneling current and substrate. The 
coupling between these states, which is mediated by the 
de-localized electrons, opens new channels for the tun- 
neling conductance which subsequently appear as steps 
in the resulting conductance measurement. 

In this Letter, we describe the STM setup in terms 
of electrons tunneling between a tip and a substrate 
onto which a sample is located. In our case the sam- 



ple comprise a spin moment which may constitute one 
or several spins. We write S = X^n^"' where the 
spin S„ is located at r„. The tunneling between 
the tip, Htip = Z^pcr^PT^crCpcr, ^ud the substrate, 

Wsub = Sk<Tcr' £k<Tdr'cj^o-Cka', occurs in presence of the 
local spin moments, and we use the model T-Lt = 
Epi.aa'n'^UTaa'irn,t)xcu.,e'^'^ + H.c.. Heref(r„,t) = 
Tg -|- T„(T ■ S„ )18l , where Tq and T„ give the tunneling 
rates for electrons not interacting and interacting, respec- 
tively, with the local spins. We include the possibility of 
magnetic tip and substrate. However, since we do not re- 
quire them to be coUinear with one another, we here have 
taken the spin quantization axis of the tip as the global 
reference frame. Due to the local nature of the spins, we 
have included a spatial dependence in the interacting tun- 
neling rate, and we use e.g. r„ = T' exp (— |r — r„|/A), 
where A is the decay length. Typically, from expansion of 
the work function for tunneling T' /Tq ~ J/U [TSj, where 
J and U are the spin-spin exchange interaction parame- 
ter and spin-independent tunneling barrier, respectively. 
For metals and semiconductors we may use T' /Tq ~ 0.1 
[To] . Here, also cr — (cr^, cr^, cr^) are the Pauli matrices. 

We note that including a Kondo-like (exchange) inter- 
action, e.g. Z^kkVcr' Cko-CTcrcr' ' S„Ck'o-', betweeu the 
local spins and the surface electrons in the substrate does 
not qualitatively alter picture for the spin moment. This 
can be understood since this interaction affects the local 
spin in a physically similar way as the interaction between 
the local spin and the tunneling current. We therefore 
omit this interaction in favor of establishing a compre- 
hensive picture of the inelastic signatures measured in 
the STM measurements. 

Using non-equilibrium technique on the Keldysh con- 
tour C, we derive the tunneling current I{t) from the 
fundamental expression 

po- 

= -|lm 5] (^P-W^-'(r-*)ck.'W)e'''"-. (1) 

pkcTcr' 

Expanding the average in the last expression according 
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to {A{t)) w (— i) J^{[A{t),HT{t')])dt' , and by converting to real times, we find that the current can be written 



pp'kk' aa' a" a 



(2) 



By expanding the tunneling operator T into its compo- 
nents To and Ti, we find that the current can be nat- 
urally written as a sum of three terms, i.e. /(t) = 
Si=o^»(^)- consider the differential conductance 

dl/dV = 'Yli^o'^hl dV of stationary source drain volt- 
ages since such conditions are predominant in experimen- 



tal situations. The first two contributions to the con- 
ductance were discussed in detail in Ref. [12] and will, 
therefore, not be further considered here. In the station- 
ary regime and at low temperatures, we can write the 
last contribution to the conductance according to 



=H 2 f ^ ^"^"^ / V^^^^^^ + eV-uj) + [l- f{eme + eV + u)) 
X {%\n{e)N(y,eF) + m(e) • M(r, e^.)]x:„(u;) 
+ [n(e)iV(r, sp) - m(e) • M(r, eF)\{Xnm{'^) + X^™(^)] 

-I- [TO^(£)A^(r,eF) -n(e)A4(r,eF)cos6'][x,7+(w) - X,|™(cj)]) ^^e. (3) 



Here, ctq — "ie^ jh is the fundamental conductance unit 
and /(e) is the Fermi function, whereas n{e) (m(e)) and 
A^(r,e) (M(r,£)) are the electronic (magnetic) densities 
in the tip and substrate, respectively, with the angle Q 
between m and M, i.e. m • M |m||M| cos^?. We have 
defined the spin-spin correlation functions e.g. Xnm('^) = 
/x-+(t)e— *di and x^™(^) = Ulm{t)e-''^'dt [20], 
where Xn^^^.t') = {-i){S- (1)8^(1')) and x^™(i,i') = 
{—i){S^{t)S^{t')) , and the bias voltage V across the 
junction. 

The conductance given in Eq. (|3]), provides a con- 
tribution to the total conductance in which signatures 
from the spin-spin correlations, or spin fluctuations, are 
present. First it is important to notice that dl2/dV is 
finite for any polarization of the tip and substrate, even 
when both are non-magnetic. Hence, under any con- 
dition, this conductance depends on fiuctuations of the 
local spin moments, here expressed by the correlation 
functions Xnmi^), xtmi^), and Xnmi^)- 

The main contribution to the differential conductance 
dl/dV is generated by the elastic tunneling processes, 
which are captured by the first contribution dIo/dV. The 
steps seen in the STM conductance experiments of local 
spins, c.f. e.g. [3 [HI [ID [H [H E] , cannot be related to 
elastic processes, because of the absence of local electron 



levels that would resonate with the bias voltage at those 
energies. Instead, the steps appearing in the dl/dV can 
be ascribed to inelastic scattering processes which open 
new channels for conductance. In the present case, such 
steps are attributed to transitions between different spin 
states. 

We may think of e.g. a local spin moment S comprising 
two coupled spins S„, n = 1,2, a spin dimer, and consider 
them to be anti-ferromagnetically coupled. The ground 
state is a spin singlet \S = 0,a = 0), while the first 
excited states constitute a spin triplet \S = Si + S2, <J — 
0, ±1). Assuming an exchange energy | J| — \Es — Et\ > 
ksT, where Et(s) denotes the triplet (singlet) energy, 
in order to prevent thermal excitations at zero bias, the 
equilibrium conductance is given by the elastic tunneling 
between the tip and the substrate only, i.e. dl/dV = 
dIo/dV. Effects from tunneling electrons scattering off 
the local spin moment averages to zero. 

The coupling to the tunneling electrons via the spin- 
spin interaction e.g. c^^cro-cr' • S^Ckcr enables, on the 
other hand, each individual spin constituting S to un- 
dergo spin-flip transitions which are assisted by spin-flips 
of the tunneling electrons. Due to this coupling, the cor- 
relation function e.g. (Tea-' ■ (S„(i)Sm(t')) • cr^'v is non- 
vanishing, in general. The spin-spin interaction, thus. 
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provides a coupling between the singlet and triplet states 
which supports transitions between them. As a result of 
these transitions, a new channel for conductance opens 
at bias voltages V > \J\/e. 

We write the Hamiltonian for the spin S„ according 
to Tin = gfJ'B^ • S„, where B is an external magnetic 
field. We account for the effective exchange interac- 
tion between the spin moments by a Heisenberg model 
Tij = —JY^n^rn ' ^his effective exchange com- 
prises a combination of e.g. direct Heisenberg exchange 
and RKKY-like (Ruderman-Kittel-Kasuya-Yosida) ex- 
change. The sign of the effective exchange parameter, 
J, may vary with distance between the spins in the 
cluster [13]. In this way we can describe clusters of 
spins with a total spin moment S. We define the re- 
sulting eigensystem {Efj, \S,a)}^^_g, for the eigenener- 
gies and eigenstates, respectively, of the model Tis = 
En nn+nj + T.n\D{Slf + E{{Slf - iSym where 
we have added the spin anisotropy fields D and E. 

The spin-spin correlation functions Xmni^)i 
and Xnmi'^) calculated in terms of the eigensystem 
for the total spin, giving 

X;F±H =(-^)27r^(z|5;F|«)(«|5±|z) 

iv 

X P{E,)[l-P{E,)]S{u + E,^Ey), (4a) 

iv 

X P{E,)[l-P{E,)]d{u + E,-E,), (4b) 

where i (v) denotes the initial (intermediate) state, 
whereas P{x) accounts for the population in the corre- 
sponding state. Eq. ([4]) provides the general qualitative 
features of the spin-spin correlation function, and shows 
that there will appear steps in dl/dV whenever the bias 
voltage matches the transition energy zt{Ei — Ey). In 
the case of a spin dimer, for instance, there appear steps 
in dl/dV when the bias voltage supports the inelastic 
transitions between the single and triplet states. Indeed, 
in the case of non-magnetic tip and substrate the present 
theory reproduces the results presented in [T7] , and shows 
an excellent agreement with the experiments in [T] [S] . 

The forms used here to describe the spin fluctuations 
given in Eq. Q correspond to the atomic limit result 
of the spin-spin correlation functions, i.e. influence from 
the tunneling current on the spin fluctuations are omit- 
ted. While this expression is a simplification, it contains 
the central mechanisms which are necessary in order to 
address the inelastic features that we are considering in 
this Letter. A more detailed description of the spin-spin 
correlation functions provides e.g. the widths (lifetimes) 
of the states involved in the transition processes, whereas 
the qualitative behavior arising from the spin-spin corre- 
lations are captured by the expressions in Eq. Q. 

Signatures of the excitations in the experimental mea- 
surements do have a finite width, which corresponds to 
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FIG. 1: (Color online) dl/dV for Fe/CuN (a) at T = 0.5 K, 
plots off-set for clarity, (fl/dV^ for (b) Fe/Pt(lll) and (c) 
Co/Pt(lll) for different widths F, at T = 4.3 K. 



that the intermediate states have finite lifetimes. Replac- 
ing the delta functions in Eq. ^ by Lorentzian funtions 
l/(a;^-|-r^), we phenomenologically include the (uniform) 
lifetime h/V for all intermediate states. In case of a single 
Fe on CuN with spin 5 = 2 [8 , and the anisotropy pa- 
rameters given in[l] we find that the Fe is weakly coupled 
(F - 10-30 ^leY) to the Cu(lOO) through the CuN layer, 
see[l](a). This value is extracted by comparing the ratio 
between the maximal and minimal conductance with the 
experimentally result (~ 1/2). Similarly, we also extract 
the widths of the single Fe (5 = 3/2) and Co (5 = 1) 
adsorbed onto Pt(lll) sur face 9;. In [1] (b), we have 
plotted d^I/dV^ for Fe (upper) and Co (lower) and we 
find the best correspondence with experiments using the 
parameters given in|T] 

In the calculations we have estimated the popula- 
tion numbers P(i?i(.y)) of the states \i{v)) involved in 
the spin-spin correlation functions, c.f Eq. Q, by 
making use of the following observation. By expand- 
ing the correlation function e.g. Xnm(^i^') terms 
of its eigenstates, we can write it as xtmi^''^') — 
E^MSn\v){v\S+\^)i-^){{d%m{dld.){t')) where 4 
(di) creates (annihilates) a particle in the state 
In the atomic limit, we can employ the decou- 
pling ((4d.)(i)(4d,)(<')) = {dlmit')){d,it)dUt')) = 
PiE.i)[l - P(£;„)]e*('^'--^")(*-*'), where the population 



TABLE 1: Anisotropy parameters, D and E, used for the 
conductance plots given in[l] the best widths, F, for the in- 
termediate states in the spin-spin correlation functions. 
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number P{Ei) — {dldi) can be estimated by using e.g. 
the Gibbs distribution P{E) = e->^'^/J2^ e"''^'- 

The rest of this Letter is devoted to inelastic spec- 
troscopy using SP-STM. It is clear from the previous 
discussion that the conductances dIo/dV and dli/dV 
only provide a change in the overall amplitude when the 
spin-polarization of the tunneling current is varied. We 
therefore focus the following discussion to the spin-spin 
correlation function, Eq. Q, in the last contribution of 
dl/dV. 

We observe in Eqs. ([S]) and Q that the transi- 
tion sequences {i\S:^\v){v\S~\i), {i\S~\v){v\S:^\i), and 
{i\S^\v){v\S^\i) are associated with different projections 
of the spin-resolved LDOS in the tip and substrate. 
The sum and difference of the first two sequences cou- 
ple to n{e)N{eF) — ni(£-) • M(£p') and mz{e)N{eF) — 
n{e)Mz{eF) cos 9, respectively, while the last sequence 
couples to n{e)N{eF) + m(e) • 'M.{£p). These different 
couplings reflect a possibility to enhance or attenuate the 
response of certain inelastic transitions by using different 
combinations of electronic and magnetic densities in the 
tip and substrate. 

In STM without spin-polarization (m, M — 0), for 
instance, the excitation spectrum can be analyzed to 
certain detail by means of applying an external mag- 
netic field e.g. along the z-direction, or the spin quan- 
tization axis of the sample, which introduces a Zee- 
man split of the levels. In order to be concrete, con- 
sider a single adatom with 5* = 1, and anisotropy fields 
D < and E = 0, which is described by the states 
\S ^ 1,(7 = 0,±1), where E^ < Eq at B = 0. In this 
system, the leading contribution to the term in dl2/dV 
containing Xnm (^) + Xmn (^) ' c.f Eq. ([3|, is proportional 
to Y.a=±i[PiEa) - P{Eo)][f{eV -E„ + Eo) - f{eV + 
Ea - Eo)], while the term containing Xnmi'^) - Xnm(^) 
vanishes. The former contribution generates steps in the 
conductance at bias voltages eV = ±{E^—Eo), separated 
by \E+i — E^i\, whenever the states |1,0) and ll,^) are 
unequally occupied. Hence, a spin split imposed by an 
external magnetic field is detectable through separate but 
equally high steps in the conductance. 

Using SP-STM opens further possibilities in the studies 
of spin systems, since then the term in dl2/dV containing 
Xnm('^) ~ Xram('^) is finite. In case of the S — 1 adatom, 
the leading contribution to this term is proportional to 
E.=±i ^[P{E,) + P{Eo) - 2P{E,)P{Eo)][f{eV ~ E„ + 
Eo) — f{e-V + Ea — Eo)]- Notice that the contributions 
to this term have opposite signs. Hence, by combining 
the electronic and magnetic densities in the tip and sub- 
strate such that m^N — nM.^ cos 6* < 0, this term leads 
to an attenuated (intensified) signal from the transition 
|l,0)(l,-hl| (|1,0)(1,-1|), while the signal is intensified 
(attenuated) when m^iV — riM^ cos 6 > 0. While the de- 
tails of the terms containing the sum and difference of the 
spin-spin correlation functions Xnm ^^'^ xtm differ from 
system to system, it is a general conclusion that the in- 



tensity of the signal from any specific transition depends 
on the magnetic densities of the tip and substrate, and 
on their relative orientation. 

We finally notice that the conductance dIo/dV van- 
ishes for nN + m • M = 0, which corresponds to having a 
half-metallic tip and substrate with their magnetic mo- 
ments in anti-parallel alignment. In this set-up also the 
conductance dli/dV — 0, which implies that the mea- 
sured signal is generated solely by dl2/dV. Moreover, 
as can be seen in Eq. this conductance only de- 

pends on the transverse components, i.e. the sum and 
difference of Xnm ^^'^ Xrtmi since the term containing 
Xnm is proportional to nN -I- m • M (= 0). Hence, de- 
spite the presence of possible thermal noise, such a set- 
up would benefit from a very low current noise since 
most of the noise would be related to the spin fluc- 
tuations that are to be measured. This can be seen 
by identifying the spin-dependent current operator with 
Si{t) = TiS{t) ■ s, where s = Epk..' 4.'^Ck.'- The 
current-current correlation function is then given by |18j 
{6i{t)Si{t')) = T^s- (S(t)S(i')) -s, where we average over 
the dynamics of the localized spins and over the ensem- 
ble of the tunneling electrons. Under the condition that 
nN -I- m • M = 0, the total dc current / is proportional 
to J^^s ■ (S(i)S(i')) ■ sd{t- t'), and since the shot 
noise is approximately (-^shot('^)) ^ the signal-to-noise 
ratio is about unity. 

In the above examples we have, for the sake of ar- 
gument, treated the local spins in the atomic limit and 
omitted mutual influences between the (electronic and 
magnetic) densities in the substrate and the adatoms. In 
reality, the adatom and the substrate affect one another 
and self-consistently creates a total effective e.g. mag- 
netic moment locally around the adatoms. The basic 
formula given in Eq. ([3| remain unchanged, however, by 
a self-consistent treatment of the system. 
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